Model calculations of the Sivers function satisfying the Burkardt 

Sum Rule 

A. Courtoy 1 , S. Scopetta 2 *, V. Vento 1 

1 Departament de Fisica Tebrica, Universitat de Valencia, 

and Institut de Fisica Corpuscular, 
Consejo Superior de Investigaciones CienUficas, 46100 Burjassot (Valencia), Spain 

2 Dipartimento di Fisica, Universitd degli Studi di Perugia, 
and INFN sezione di Perugia, via A. Pascoli 06100 Perugia, Italy 

Abstract 

It is shown that, at variance with previous analyses, the MIT bag model can explain the available 
data of the Sivers function and satisfies the Burkardt Sum Rule to a few percent accuracy. The 
agreement is similar to the one recently found in the constituent quark model. Therefore, these 
two model calculations of the Sivers function are in agreement with the present experimental and 
theoretical wisdom. 
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The study of the partonic properties of transversely polarized hadrons will answer crucial 
questions on their structure, such as their relativistic nature and their angular momentum 
content. Experiments are progressing fast, motivating a strong theoretical activity [1]. One 
of the quantities under scrutiny is the Sivers function, the object of this study. Semi-inclusive 
deep inelastic scattering (SIDIS), i.e. the process A(e, e'h)X, with the detection in the final 
state of a hadron h with the scattered electron e', is one of the processes to access the 
transversity parton distributions (PDs). For several years it has been known that SIDIS off 
a transversely polarized target shows the so called "single spin asymmetries" (SSAs). It can 
be shown (see, i.e., Ref. and references therein), that one of the mechanisms generating 
the SSAs is governed by the Sivers function jsj]. The latter describes the amplitude of 
modulation of the number density of unpolarized quarks in a transversely polarized target 
due to the correlation between the transverse spin of the target and the intrinsic transverse 
parton momentum. The Sivers function is a Transverse Momentum Dependent (TMD) PD, 
denoted kx), where x is the Bjorken variable and kx is the transverse momentum of 

the parton Q. It is a time reversal odd object [jj and for this reason, for several years, it was 
believed to vanish. However, this argument was invalidated by a calculation in a spectator 
model, following the discovery of Final State Interactions (FSI) at leading- twist, i.e., not 
kinematically suppressed in DIS [4]. The current wisdom is that a non- vanishing Sivers 
function is generated by FSL technically represented by the gauge link in the definition of 
TMD parton distributions [5[. Recently, the first data of SIDIS off transversely polarized 
targets have shown a strong flavor dependence of the Sivers mechanism 6J. Complementary 
experiments on transversely polarized 3 He target, addressed in [3], are being performed at 
JLab [8]. Parameterizations of f^(x, kx) are available 0, Q, , and new data are expected 



soon. From the theoretical point of view, a mode 



independent constraint on calculations 



of f 1T ^(x,kx) is the Burkardt Sum Rule (SR) 12j. It states that the average transverse 
momentum of all the partons in a hadron, (kx), which can be defined through f^(x, kx), 
has to vanish. If the proton is polarized in the positive x direction, the Burkardt SR reads: 

E <*?> = ~ fdx f dkx^f^ix, k T ) = 0. (1) 

Q=u,d,s,g.. ^° J 

Given the present situation of increasing experimental activity, estimates of f^(x, kx), 
subject to solid theoretical constraints, can be very useful. Since a direct calculation in 
QCD is not yet feasable, this quantity has been calculated in several models: a quark- 



diquark model 4], Ll3|] ; the MIT bag model, in its simplest version 



4| and introducing 



an instanton contribution [15j ; the Constituent Quark Model (CQM) [16]. To distinguish 
between the model estimates, data and model independent relations, such as the Burkardt 
SR, can be used. In all the models used so far the total momentum of the proton is carried 
by the quarks of flavor u and d. According to Eq. ([T]), this implies that the magnitude of 
f^j? for Q = u and d has to be similar and the sign has to be opposite. This is also the 



trend of the parameterizations of the data 



In different versions of the diquark model 



the magnitude of the d contribution is much smaller than that of the u. In the MIT bag 



model [14j u and <i-quark contributions of opposite sign are found to be proportional by a 



factor of -4. Even in the modified version of the MIT bag model of Ref. 151 ]. the Burkardt 
SR is not fulfilled. On the contrary, in the CQM, we found a satisfactory description of 



the data and therefore the calculation fulfills the Burkardt SR at the 2% level [16]. This 
puzzling situation deserves to be investigated. To this end, we next analyze the MIT bag 
model calculation to understand the origin of the discrepancy with the CQM calculation. 
One should realize that, in the CQM, even if a pure S-wave description of the proton is 
used, i.e., a pure SU(6) wave function, we are able to reproduce the gross features of the 
data. The same SU(Q) spin-flavor structure is used in the bag calculation of Refs. [l4, 2] 
and no agreement with the data is found. This situation is in contradiction with previous 
calculations of other PDs in the bag model [3] and in the CQM [3] which both have been 
able to reproduce the gross features of the data. 

If the proton is polarized in the positive x direction, the Sivers function can be written, 
in a helicity basis for the proton, as la, Il7 ] 



f^f(x,k T ) 



29? 



M 
4fc„ 



dj-d 2 j T 
(2vr) 3 ' 



-i(x^-P+-^ T -k T ) 



(ps z = i\o q \ps z 



(2) 



where Q = ipa(0Ji , &)£t (oo, f )7 + £ (oo, O)Vg(O) , V>e(0 is the q uark field^jC^oo, £~ 
is the gauge link 
to calculate Eq. 



11] . In the following, the framework and the notation of Ref. 14j are used 
2]) in the MIT bag model Q. By expanding the gauge link to next to 



eading order, inserting in Eq. (T5]) the bag model wave function in momentum space, <p(k) 
2o| , and using for the definition of the quark helicity and momentum labels the ones in 
Fig. [H /{i^ s (x, k T ) can be written as 
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k,m 2 



k — q, m 1 



ki,m,4 (O k 3 ,m 3 




FIG. 1: The contributions to the Sivers function in the present approach. The graph has been 
drawn using JaxoDraw [24]. 
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This equation corresponds to Eq. (17) in Ref. [14J] modified to follow the Trento conven- 
tion [17|] which implies an additional factor of 1/2. Here g is the strong coupling constant, 
Cf4 = T«T^(PS Z = l\4 mi ^ Qm AArnJPS z = -1) with {m} = mi ,m 2 ,m 3 ,m 4 ; M is 
the proton mass, E p its energy, 6g m is the annihilation operator for a quark with flavor Q, 
helicity m, and color index i, and is a Gell-Mann matrix. In turn, the k 3 integral can be 
written as 

d 3 k 



v3 ,„t 



(27T) 



3^3(^)7 7>m*(fcs-?l) = F ma (q ± ) 5 m3m4 +H m3 (q ± )5, 



m 3 ,-m4 ; 



(4) 



with 



^m 3 (9-1 



c 

71 



tptg + k^t^Q jk 3 + A; 3 tg/fc 

+ {h-k' 3 + iv Z d m3 )t\tl/{k' 3 h) , 

ff m3 (gl) = ^=J d% [{ikl - k%d m3 ) t'X/h - (ik'i - k'i d m3 ) t 3 t 3 /k> 3 

+ (vyd m3 + tv x )tlt 3 /(hk' 3 )\ , 



„^3/3 



(5) 



(6) 



where k 3 



I ^3 1 5 ^3 — I ^3 1 J ^Ti3 — (^m 3 i 



"13,- 



1 , v 



Q±><k 3 , k' 3 = k 3 - q, C 



1Quj a / {ii 2 il{uj){ijj — l)Mp), with a; being the bag model mode [20| and the function t 
U(k 3 ), tf = U(hQ, i = 0, 1 are defined in Q- 



In Ref. [lij], only the first term in the r.h.s. of Eq. (j3J) is calculated. This term corresponds 
to the helicity conserving contribution (1713 = 1714) associated with the quark in the lower 
part of Fig. [TJ This result only applies if the integral is performed taking q± along the z 
direction. However, this is incorrect in the present case. As in any DIS process, the direction 
of the virtual photon determines the operator structure, i.e. 7 + in here. The fixing of the 
photon direction leads to a 73 matrix in this operator. Therefore we do not have anymore 
the freedom to choose z as the direction of the exchanged gluon, which must lie in the (x, y) 
plane. Besides, one can check that the integral Eq. (J4]) does depend on the direction of q±. 
Moreover, if the findings of Ref. [jjl were correct, it would mean that a helicity flip could 
occur only for the quark which interacts with the photon and not for the quark in the lower 
part of Fig. [U a restriction which does not have any physical motivation. Thus the present 
calculation differs from the previous one in that we take into consideration both terms of 
Eq. By the same argument, the expression tp\ ril {k — cf_j_)7 7 + y? m2 (fc) in Eq. ([3]) also 



contains both helicity-flip and non-flip terms. One should notice that, in Ref. [15l |. a helicity 
flip term (with = —1714) has been found to contribute to the Sivers function. In that 

)aper, instanton effects have been added to the pure MIT bag model calculation of Ref. 

14J and the presence of this helicity flip term is due solely to the instanton contribution. 
However, in the calculation of Ref. [l5(, the term depending on S m3 - m4 in Eq. (4), associated 



to the perturbative one-gluon exchange, should appear and has not been considered 



2l|. 



It is interesting to realize that, in a completely different scenario, the CQM calculation 
satisfying the Burkardt SR of Ref. [if]], a contribution is found for helicity conserving and 
helicity flip of the quark in the lower part of Fig 1, as it happens in our MIT bag model 
calculation. 

Evaluating the matrix elements for the valence quarks and assuming an SU(6) proton 
state in Eq. ([3]), one gets 

k ± ) = -t ML c> J ^y 2 ^[C Q + Y(q ± , k T ) + C+-U(q ±J hr)] , (7) 
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FIG. 2: The quantity f^p q (x). Eg. (|10p . for the u and d flavour. The dashed curves are the 
results of the approach of Ref. [14], the full ones those obtained here. 



with 



Y(q ± , kr) = [k' y t%/k' - kyht'Jk - v y t 1 t[/ (kk')] 



d 3 k. 



t^t^ + k^tlt^ /k 3 + k^t^t^/k^ 



^3 
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3+'3 ni 



,x +3j.'3 



(9) 



where k = \k\, k' = \k'\, U = t^k), t\ = t^k'), i = 0,1 and C~+ = -16/9 {C^ + = 4/9), 
C+ _ = -4/9 (CJ+ = -8/9) for Q = u{d). Let's recall that in Ref. Q only the first term of 
Eq. ([7]), proportional to contributes to /fy 2 . It is therefore found that = — • 

Notice that, in order to calculate /f^ 2 , one is using a two-body operator associated with FSI 
and therefore one should not expect a proportionality between the u and d results. On the 
contrary, in the calculation of conventional PDs, in any SU(6) model calculation, the used 
operators are of one-body type and therefore the results turn out to be proportional [18 ]. 
Numerical results are shown in Figs. [2] and [3] for the first moment of /f^ 2 , i.e. 



X(l)< 



d 2 k T 3^fiT Q ^M) 



(10) 



In Fig. [2] the dashed curves are the ones obtained in Ref. (cf. Fig. 4 in Ref [l^] adapted 
to the Trento convention [l3], i.e., reduced by a factor of two). The obtained value for the 
Burkardt SR Eq. ([T]) turns out to be 8.73 MeV. To have an estimate of the quality of the 
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FIG. 3: The same as in Fig. 2, after NLO evolution (see text). The patterned area represents the 



1 — cr range of the best fit of the HERMES data proposed in Ref. 



Ill- 



agreement of this result with the SR, we consider the ratio r = ({k%) + (k%))/({k%) — (k%)) , 
obtaining r ~ 0.60, i.e., the Burkardt SR seems to be violated by 60 %. The full curve in 
Fig.[5]is the result of the present calculation. Clearly, the d contribution becomes comparable 
in magnitude to the u one. The obtained value for the Burkardt SR is -0.78 MeV and 
r ~ 0.05, i.e., it is only violated by 5 percent. These results are comparable in quality 
to those obtained for the CQM 161 ]. restoring the approximate agreement between the 
two schemes. In order to compare the results with the data, one should perform a QCD 
evolution from the experimental scale, which is, for example for the HERMES data, Q 2 = 2.5 



GeV 2 



221 ] . Unfortunately, the evolution of TMDs is still to be understood, although recent 



developements can be found in Ref. 231 ] . In order to have an indication of the effect of 
the evolution, we evolve at NLO the model results assuming, for the moments of the Sivers 
function, Eq. ffTUj) . the same anomalous dimensions of the unpolarized PDFs, as we did 
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FIG. 4: The same as in Fig. 3, but comparing with the parameterization of the data proposed in 
Ref. [9] (patterned area). 

in Ref. [if]] for the CQM calculation. The parameters of the evolution have been fixed in 
order to have a fraction ~ 0.55 of the momentum carried by the valence quarks at 0.34 
GeV 2 , as in typical parameterizations of PDs , starting from a scale of /Xq ~ 0.1 GeV 2 
with only valence quarks. The results in Fig. [3] and H] show an impressive improvement of 
the agreement with data once the full contribution of Eq. ([7j) is taken into account. The 
data are described rather well for both flavors. Comparing this encouraging outcome with 
that of Ref. [16J, one can notice that the Burkardt SR is better fulfilled in the CQM. Most 
probably this has to do with the fact that the Burkardt SR is associated with transverse 
momentum conservation and, in the MIT bag model, the proton wave function is not an exact 
momentum eigenstate. In closing, we can say that, for the first time, it has been established 
that correct model calculations provide phenomenological successful interpretations of the 
Sivers function, which are consistent with each other. 
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